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I.  INTRODUCTION 

This  report  describes  a digital  simulation  for  the  evaluation 
of  chaff  decoy  placement  and  related  antiship  missile  defense  problems. 

The  simulation  is  the  latest  in  a series  of  programs  developed  by  Systems 
Consultants,  Inc.  (SCI)  for  the  U.S.  Navy;  this  one  was  developed  for  the 
Naval  Surface  Weapons 'Center  White  Oak  Laboratories  under  contract 
N60921-77-C-0089.  The  purpose  of  that  contract  was  to  investigate  logistic 
requirements  for  the  MK  36  and  MK  182  Super  Rapid  Bloom  Off board  Chaff 
(SRBOC)  decoy  system. 


This  document  represent 


)inion  of  Systems  Consultants, 


Inc.  The  views  and  conclusions  set  north  herein  are  those  of  the  author 
and  should  not  be  interpreted  as  necessarily  representing  the  official 
policies,  either  expressed  or  implied,  of  the  Naval  Surface  Weapons  Center 
or  the  U.S.  Government. 
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II.  OVERVIEW  OF  THE  SIMULATION 


A.  BACKGROUND 

A decoy  is  a device  which  is  designed  to  present  an  appearance 
which  is  indistinguishable  from  that  of  a true  target  when  viewed  by  a 
specified  sensor,  so  that  deception  of  the  sensor  or  its  operator  can  be 
accomplished.  The  degree  to  which  the  decoy  must  resemble  the  true  target, 
in  terms  of  all  of  its  measurable  characteristics,  depends  upon  the 
attributes  of  the  sensor.  For  example,  a sensor  which  measures  both  the 
signature  and  motion  characteristics  of  its  target  cannot  be  reliably 
deceived  by  a decoy  which  differs  significantly  from  a true  target  in 
either  respect. 

Decoys  have  been  in  use  since  the  earliest  days  of  radar.  The 
ways  in  which  decoys  can  be  used  to  defeat  weapons  correspond  to  the  sequence 
of  operations  of  weapon  targeting.  During  World  War  II,  for  example, 
chaff  decoys  were  used  by  the  allies  against  German  air  defenses  in  two 
ways:  corridor  chaff  was  employed  to  mask  entire  flights  of  bombers  to 

deny  information  on  the  number  of  aircraft  involved;  chaff  was  also  used 
to  degrade  the  tracking  capabilities  of  the  German  gun-laying  radars. 
Currently,  decoys  are  being  developed  to  counter  anti  ship  missiles  (ASMs). 
Against  these,  there  can  be  as  many  as  three  modes  of  decoy  operation: 

(1)  Target  Selection  Confusion  of  the  Operator  (TSCO),  wherein 
the  purpose  of  the  decoy  is  to  cause  the  weapon  to  be 
launched  against  a decoy  rather  than  a true  target  or  to 
deny  targeting  data. 

(2)  Initial  Acquisition  Decoying.  Many  guided  weapons, 
particularly  long  range  ones,  do  not  track  their  targets 
continuously  from  launch,  but  fly  under  the  control  of 
an  autopilot  for  some  specified  period,  after  which  the 
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target  must  be  acquired.  The  purpose  of  initial 
acquisition  (lA)  decoying  is  to  cause  the  missile 
to  acquire  a decoy  rather  than  a true  target. 

(3)  Break  Lock  Decoying.  Once  a missile  has  acquired  a 
true  target,  break  lock  (BL)  decoys  can  be  used  to 
degrade  the  tracking  accuracy  of  the  missile. 

The  foregoing  example  thus  represents  instances  of  TSCO  and  BL  decoying. 

The  efficacy  of  a decoy  depends  upon  satfsfaction  of  two 
criteria:  indistinguishability  from  a true  target  and  placement  relative 

to  the  true  target  and  the  sensor.  The  decoy  must  be  indistinguishable 
from  a true  target  in  the  sense  that  the  victim  sensor  cannot  reject  the 
decoy  on  the  basis  of  received  data.  The  decoy's  position  is  important 
because  sensors  accept  only  those  signals  which  emanate  from  the  volume 
of  space  where  true  targets  are  expected  to  be.  Sensor  characteristics 
may  differ  according  to  their  mode  of  operation  (pre-launch  targeting, 
acquisition  of  the  target  by  the  weapon,  and  terminal  homing)  and  so  the 
requirements  for  decoy  characteristics  and  placement  may  also  vary  accord 
Ing  to  decoy  mode  (TSCO,  lA  or  BL).  Since  both  the  indistinguishability 
and  placement  criteria  must  be  satisfied  for  the  decoy  to  be  effective, 
they  can  be  evaluated  independently  in  most  instances.  This  principle  is 
central  to  the  utility  of  the  digital  simulation  described  in  this  report 
which  is  used  to  evaluate  decoy  placement  in  the  BL  mode  against  radar 
guided  ASM. 

The  decoy  placement  criteria  can  always  be  satisfied  if  there 
Is  sufficient  selectability  in  decoy  placement  relative  to  the  defended 
ship  (true  target)  and  the  victim,  and  if  the  timing  of  the  decoy  launch 
Is  appropriate.  However,  practical  limitations  to  decoy  placement  will 
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always  exist.  Some  of  these  limitations  arise  as  a result  of  uncertainity 
as  to  the  location  (range,  altitude  and  direction),  speed,  direction  of 
motion  and  detailed  characteristics  of  the  ASM.  Others  result  from  the 
design  limitations  of  the  decoy  system:  in  order  to  minimize  equipment 

cost  and  size,  shipboard  decoy  launchers  employ  fixed  train  and  elevation 
angles;  the  decoys  themselves  are  ballistic  and  use  a single  propulsion 
system.  Consequently,  the  design  of  decoy  systems  to  achieve  maximum 
placement  effectiveness  and  the  evaluation  of  the  placement  effectiveness 
of  a decoy  system  are  relatively  complicated  problems.  The  simulation 
described  subsequently  in  this  report  was  developed  specifically  to 
address  these  problems. 

B.  SCOPE 

Decoy  effectiveness  against  a particular  ASM  depends  upon  the 

decoy  placement  geometry  and  also  upon  the  performance  characteristics  of 
the  decoy  (e.g,,  its  radar  cross  section).  Consequently,  it  is  possible 


to  determine  decoy  effectiveness  only  through  consideration  of  both 
performance  and  placement.  Since  these  are  independent  of  each  other  for 
most  applications  (and  in  particular  for  the  threat  model  used  herein) 


it  is  possible  to  write  the  following  expression  for  decoy  effectiveness: 


^ ■ ^placement  ^performance 

where  ^placement  the  placement  effectiveness  and  ^performance 
is  the  performance  effectiveness,  i.e.  the  probability  that  the  missile 
will  home  on  the  decoy  given  that  the  decoy  is  placed  within  its  field 
of  view.  The  decoy  placement  simulations  have  been  designed  to  calculate 


^placement,  and  therefore  yield  r],  the  decoy  effectiveness,  only  under 
the  condition  that  ^performance  equal  to  unity.  That  is,  the  simulation 
gives  the  upper  bound  on  decoy  effectiveness. 
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The  results  from  the  simulation  are  further  limited  because 
of  designed-in  assumptions  about  the  nature  of  the  engagement.  Placement 
effectiveness  depends,  inter-alia,  upon  the  range  and  direction  of  the 
threat  at  the  time  of  decoy  launch,  the  speed  of  the  ship,  and  the  speed 
and  direction  of  the  true  wind.  The  simulation  yields  average  placement 
effectiveness  values  for  a given  ship  speed,  decoy  launcher,  and  decoy 
against  a specified  threat  by  sampling  from  a three  dimensional  space 
which  is  bounded  by  minimum  and  maximum  threat  directions,  wind  directions, 
and  wind  speeds.  The  first  of  these  is  uniformly  distributed,  i.e.  any 
value  between  the  bounds  is  equally  likely.  The  wind  conditions  (true 
speed  and  true  direction)  can  be  chosen  in  either  of  two  ways.  The 
preferred  method  is  to  choose  wind  direction  from  a uniform  population, 
and  to  select  the  speed  from  a parent  population  selected  to  closely 
match  empirical  data  (Long  et.  al.,  1965).  The  second  method  involves 
successively  incrementing  true  wind  speed  or  true  wind  direction  on  a case- 
by-case  basis,  so  as  to  arrive  at  a sample  which  includes  a predetermined 
matrix  of  true  wind  speed  and  true  wind  direction  values.  For  obvious 
reasons,  this  option  is  referred  to  as  determistic.  As  a result  of  the 
sampling  method,  the  average  results  of  the  simulation  cannot  be  applied 
to  any  one  scenario,  but  represent  the  placement  effectiveness  achieved 
over  a large  number  of  scenarios  which  occur  under  varying  conditions. 

Since  the  simulation  has  provisions  for  detailed  reports  giving  the  place- 
ment results  on  a sample  by  sample  basis,  it  can  also  be  used  in 
evaluating  decoy  placement  performance  in  any  particular  scenario. 

C.  STRUCTURE 

The  simulation  represents  the  embodiment  of  mathematical  models 


and  evaluation  criteria  in  software.  There  are  three  elements  to  the 
mathematical  model:  the  decoy,  ship  and  threat.  The  mathematical  model 
describes  the  salient  features  of  each  element  and  their  motion  relative 


to  each  other.  The  evaluation  criteria  consist  of  sets  of  precise 
conditional  statements,  the  outcomes  of  which  are  classified  as  pass  or 
fail  with  respect  to  various  success  determinants.  Implemented  in  software, 
the  mathematical  model  and  evaluation  criteria  are  used  with  statistically 
independent  samples  of  the  appropriate  variables  to  perform  a time  line 
analysis  of  an  engagement.  By  manifold  replication  of  this  process, 
average  placement  effectiveness  values  are  obtained  (Monte  Carlo  method). 

The  software  also  includes  report  generators  of  two  types:  summary  and 
detailed.  The  summary  reports  yield  aggregated  results  (e.g.  average 
placement  effectiveness)  while  the  optional  detailed  reports  yield  results 
(values  of  the  random  variables  and  pass/fail  indication)  on  a case  by 
case  basis. 

The  mathematical  model  upon  which  the  simulation  is  based  is 
discussed  in  Section  III  of  this  report.  Section  IV  describes  the 
evaluation  criteria,  and  Section  V includes  FORTRAN  IV  software  and 
interactive  user's  guide  for  the  simulation  program. 


III.  THE  MODEL 


The  simulation  includes  sets  of  equations  for  the  coordinates 
of  the  ship,  decoy  and  threat  in  two  reference  frames:  the  earth  frame 
(so  called  because  it  is  at  rest  on  the  surface  of  the  earth)  and  the 
missile  frame  (which  is  at  rest  with  respect  to  the  missile).  These 
reference  frames  are  connected  by  means  of  a time  dependent  coordinate 
transformation:  at  any  instant  the  transformation  consists  of  a trans- 

lation plus  a rotation  of  coordinates.  Once  the  transformation  is  effected, 
the  coordinates  of  any  point  are  given  in  both  frames.  By  referencing  the 
missile  frame  to  its  line  of  flight  direction,  the  angles  from  boresight 
to  any  r n space  can  be  determined,  that  is,  it  is  deteminable 
wheth  nt  lies  within  the  missile  seeker  antenna  beamwidth. 

Natu  j,  it  is  also  easy  to  determine  the  angle  subtended  by  an  object, 
and  also  distances.  In  the  following  subsections  the  equations  of  motion 
for  the  ship,  decoy  and  missile  are  developed.  Field  of  view  and  range 
considerations  are  discussed  in  section  IV,  as  they  are  naturally  related 
to  the  evaluation  criteria. 

A.  SHIP 

For  the  purposes  of  these  simulations  the  ship  is  represented 
as  a line  of  length  2BSHD  (^ow-£tern  Half  Distance).  Its  motion  is  uniform 
in  the  earth  frame,  that  is,  it  travels  at  a constant  speed  in  a 
constant  direction.  The  reference 
direction  is  chosen  to  be  along  the 
x-axis  of  the  earth  frame;  at  time 
T*0  the  midpoint  of  the  ship  is 
located  at  the  origin  of  the  co- 
ordinate system,  as  shown  in  Figure 
(1).  (In  the  figure,  the  vertical 


FIGURE  1 


direction  is  normal  to  the  page).  The  equations  of  motion  for  the  ship's 
midpoint  are  thus: 


Xs=VsT 


(1) 


Ys=Zs=0 

while  the  coordinates  of  the  bow  and  stern  are  shifted  along  the  x-axis 
by  +BSHD  and  -BSHD,  respectively; 


hs‘h  * 85H0  (Vss=ZSB=0) 

Xss=X$  - BSHD  (YsS=ZsS=0) 

It  has  been  determined  that  the  ship  can  be  accuately  modeled 
in  terms  of  a line  rather  than  a three  dimensional  body  for  the  following 
reasons: 


(1)  The  missile  track  point  is  generally  within  a few  meters 
of  the  ocean  surface,  so  that  it  is  not  necessary  to  re- 
present the  ship  in  terms  of  a vertical  extent. 

(2)  The  ship  is  very  long  in  comparison  to  its  width 
(typical  destroyers  have  L/B  ratios  of  about  9),  so  as 
a consequence  the  width  of  a ship  contributes  little 
to  its  angular  extent  when  viewed  from  distances  of 
practical  interest.  Further,  at  these  distances  the 
angle  subtended  by  the  ship's  width  (when  viewed  bow- 
on  or  from  astern)  is  typically  very  much  smaller  than 
the  beamwidth  of  the  missile.  Thus  it  is  not  necessary 
to  complicate  the  model  by  treating  the  ship  as  even 

a two  dimensional  body. 

It  is  worthwhile  to  comment  that  the  modeling  of  the  ship  even 
in  terms  of  such  a simple  entity  as  a line  has  very  little  effect  upon  the 
outcome  of  the  simulation  analysis  except  in  the  case  where  a very  high 
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performance  missile  (narrow  beamwidth  and  range  gate)  is  assumed;  this 
applies  likewise  to  the  decoy.  The  initial  simulation  which  was  used  by 
SCI  to  evalaute  decoy  placement  effectiveness  represented  both  the  ship 
and  the  decoy  as  points,  yet  exhibits  the  same  sort  of  qualitative  behavior 
as  does  this  one  with  respect  to  the  more  fundamental  of  the  success 
criteria. 

B.  DECOY 

A decoy  is  treated  within  the  model  as  an  object  which  in  effect 
comes  into  existence  at  some  precalculated  time  and  position  relative  to 
the  ship.  This  reflects  the  need  to  model  the  times  of  flight  and  cross 
section  development  of  the  decoy,  (These  times  are  called  Tf  and  Tb). 
Suppose  that  the  decoy  is  launched  at  time  T=0  from  a launcher  located 
a distance  d[_  forward  of  amidships  and  pointed  at  an  angle  4>  off  the  bow, 
as  shown  in  Figure  2.  (Note  that  angles  are  measured  according  to  the 
mathematical,  not  nautical,  convention.)  Let  be  the  horizontal  distance 
of  decoy  deployment  and  ho  be  the  ^ ^ 

altitude  of  decoy  deployment.  Since 

the  decoy  initially  has  a velocity  com-  , 

ponent  Vs  in  the  X direction,  its  position 

. A 

at  time  T=Tf  will  be  given  by  ^ ^ ^ 

-)l  k- 

Xc“dL+RcCOS({)+VsTf  FIGURE  2 

Yc*Rcsin4)  (3) 


FIGURE  2 


Since  the  decoy  is  airborne,  it  will  be  subject  to  vertical  motion  under 
the  influence  of  gravity  and  horizontal  motion  under  the  influence  of  the 
true  wind,  as  measured  in  the  earth  frame.  Let  6 be  the  angle  of  the 
true  wind  as  measured  off  the  x-axis  in  Figure  2 and  let  the  true  wind 
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speed  by  Vyj.  Suppose  also  that  the  decoy  quickly  reaches  its  terminal 
velocity  of  fall  Vcy.  Then  the  coordinates  of  the  decoy  at  time  T>Tf  are 
XQ=dL+RQ  cos(()  +VjT^  -Vw  cos0 

Yc=Rcsin(()-VyTsine  (4) 

Zc=ho  -Vcv(T-Tf) 

The  decoy  is  treated  as  a sphere  of  radius  Dec  centered  aL  the  coordinates 
given  by  (4). 

The  horizontal  motion  of  the  decoy  relative  to  the  ship  will  be 
determined  by  the  relative  wind  (wind  speed  and  direction  as  measured 
by  the  ship  - See  Schiff,  et.  al . , 1975).  Since  the  relative  wind 
conditions  are  of  tactical  importance,  the  calculation  of  relative  wind 
speed  and  direction  are  included  in  the  simulation.  The  relative  wind 
speed  is  given  by 


while  the  direction  is 

RWA=sin-l  (VyjSin  (6)/Vr|^) 

This  second  equation  has  two  admissible  roots  which  differ  by  180°  - 2RWA. 
The  physically  correct  root  is  determined  by  considering  which  quadrant 
the  true  wind  lies  in  and  its  speed  relative  to  that  of  the  ship.  The 
algorithm  for  accomplishing  this  is  discussed  in  Appendix  A. 

C.  MISSILE 

The  missile  model  utilized  in  the  simulation  has  been  designed 
for  maximum  generality,  to  minimize  reliance  upon  the  need  for  detailed 
intelligence  information.  The  general  characteristics  of  the  model  are: 

(1)  Angle  gating  so  as  to  process  information  emanating  within 
a cone  centered  on  the  missile  boresight,  of  central  angle 
2 BWHA  (BeamWidth  Half  A^ngle). 

(2)  Range  gating  to  further  restrict  the  volume  of  space  from 
which  processed  signals  can  emanate.  The  width  of  this 
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gate  is  2RGHW  (Range  £ate  H^alf  Width).  The  range 
gate,  like  the  angle  gate,  is  assumed  to  be  geometrically 
centered  on  the  traget. 

(3)  Fixed  ground  speed  Vj 

(4)  Fixed  altitude  Zj 

(5)  Proportional  (Pursuit)  navigation  in  the  azimuth 
plane  (but  see  below!). 

For  reasons  of  economy  of  simulation  design  and  execution, 
the  flight  of  the  missile  is  decomposed  into  two  segments:  prior  to  and 

subsequent  to  detection  of  the  decoy.  From  its  initial  point  the 
missile  moves  along  a straight  line  toward  the  estimated  point  of  impact 
with  the  ship  midpoint  until  the  decoy  is  detected.  Upon  satisfaction 
of  the  decoy  detection  criterion,  the  missile  turns  toward  the  decoy  at 
a rate  which  is  proportional  to  the  azimuth  angle  off  boresight  to  the 
decoy. 

The  two  segments  of  the  missile  dynamic  model  represent 
different  navigational  modes:  lead  and  proportional.  This  is  in 

contrast  to  real  missiles,  which  would  employ  a single  navigational 
mode.  However,  the  differences  between  the  flight  paths  for  the  two 
modes  depend  in  the  main  upon  the  relationship  between  the  target  and 
missile  velocities.  Since  the  missile  speed  is  at  least  an  order  of 
magnitude  greater  than  those  of  either  the  ship  or  the  decoy,  the 
difference  between  the  flight  profiles  for  the  two  modes  is  insignificant, 
and  thus  the  simplification  made  in  the  model  is  entirely  justified. 

At  time  T=0  the  missile  is  a distance  Rq  from  the  ship  at 
bearing  y.  The  missile  is  moving  in  the  xy  (horizontal)  plane  with  speed 
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Vf,  while  the  ship  moves  along  the  x-axis  of  the  earth  frame  with  speed 


V5.  Let  Tg  be  the  time  of  impact;  the 
missile  will  move  a distance  V-rTg  and  the 
ship  will  move  a distance  V^Tg  in  the  time 
interval  (0,Te).  Using  the  law  of  sines 
B=sin  {Vs  siny}  (5) 

W 

subject  to  Vs<V-[-,  and 

^ _ Ro  siny  (6) 

®”  VjSin  (B+y) 

With  these  relationships  the  motion  of  the 
given  by 


''S'e 

FIGURE  3 

missile  in  the  xy  plane  is 


Xy=Rg  cosy  (l-T/Tg)+VsT 
Yy=Ro  siny  (1-T/Tg) 
for  times  T<Tg. 

The  missile  frame  of  reference  is  defined  as  that  which  is  at 
rest  with  respect  to  the  missile,  with  its  origin  corresponding  to  the 
missile  location  and  its  x-axis  directed  along  the  missile  line  of  flight, 
as  shown  in  Figure  3.  The  transformation  between  the  earth  frame  and  the 
missile  frame  is  accomplished  by  first  making  a translation 


X'  = X-Xy 
r = Y-Yy 

followed  by  a rotation  in  a counterclockwise  direction 
X"  = -X'  cos  (6+y)-Y'sin  (B+y) 

Y"  = X'sin  (B+y)  -Y'cos  (B+y) 


X"  « (Xy-X)cos(B+y)+(YT-V)sin  (B'h')  (7) 

Y"  - (X-Xy)sin(B+y)+(Yy-Y)cos  (B+y) 


The  coordinate  transformation  also  includes  a translation  in  the  vertical 


direction: 

Z"*Z-Zt  (8) 

During  the  decoy  tracking  segment  of  the  missile  flight,  the 
missile  trajectory  in  the  horizontal  plane  is  determined  by  the  position 
of  the  target  relative  to  the  missile  line  of  flight.  In  particular,  the 
direction  and  rate  of  turn  in  the  azimuth  plane  are  related  to  the  angle 
off  bores ight  by 

da  *a=-Ka  (9) 

3t 

where  a is  the  rate  of  turn,  a is  the  azimuth  angle  off  boresight, 
and  K is  the  navigational  system  constant.  Equation  (9)  describes  a 
simple  proportional  navigation  system.  The  constant  k is  related  to  the 
acceleration  limit  of  the  airframe^  and  in  turn  to  the  sensor  beamwidth. 

Suppose  at  some  instant  tQ  the  missile  is  flying  at  speed 
in  the  positive  x-direction  in  some  reference  frame;  at  that  time  an 
azimuth  error  angle  a exists.  According  to  the  navigation  system 


equation  (9),  the  missile  line  of  flight  is  rotated  thru  the  angle  ic  a 6t 
during  the  infinitesimal  time  interval  6t,  so  that  at  time  t = tg  + 6t 
the  X and  y velocity  components  of  the  missile  are 


Vy  (to+6t)  = Vt  cos  Ko6t  Y 


and 


Ay  = Lim  Vy  (to+6t)  - Vy  (to) 

6t-o  6t 

= Lim  Vt  sin  <a6t  = KaVt 
6t-K)  6t 

Consequently,  if  %ax  maximum  azimuth  angle  error  which  can  be 

accepted  and  amax  is  the  airframe  maneuver  limit, 

®max  “ ^ “max  '^t 

so  that 

K = ^max  (10) 

“max  ^t 

In  order  to  calculate  k,  is  generally  taken  to  be  approximately 
equal  to  the  scanned  half  bearawidth  limits  of  the  missile  radar  antenna. 

On  the  basis  of  some  further  simplifications,  the  required  maneuver 
limit  can  be  calculated  based  upon  the  missile  field  of  view  geometry 
so  that  < can  be  uniquely  determined  on  the  basis  of  the  missile  speed, 
search  beamwidth  and  minimum  target  acquisition  range  (Schiff,  1977). 

The  algorithm  for  the  tracking  of  the  decoy  (subsequent  to  the 
break  lock  decision  point)  is  as  follows.  At  the  moment  of  lock  transfer, 
the  missile  coordinates  in  the 
earth  frame  are  (Xt,Y^);  the 
ship  is  located  at  (Xs,0)  and 
chaff  cloud  coordinates  are 
(Xc»Yc)*  At  that  instant,  the 
missile  is  proceeding  toward 
the  ship  as  shown  in  Figure 
5.  The  error  angle  c^is  deter 
mined  by  first  determining  the 
orientation  of  the  missile  coordinate  frame. 


the 


FIGURE  5 


y'  * tan  (Yt/Xt) 


e'  = sin  (XsCosyV  (Xt-Xs)^+Y|)  ) 

The  angle  of  rotation  between  the  earth  frame  and  the  missile  frame,  for 
the  missile  in  the  first  or  second  quadrants,  is  180°  - (3'+y‘)- 
Accomplishing  a translation  of  coordinates  to  (X^,  Yt)  followed  by  this 
rotation,  the  coordinates  of  the  decoy  are  calculated.  (Cf  equation  (7).) 
The  error  angle  ap  is  then  given  by:  ap  = tan  (Y^/X^) 

In  the  example  of  Figure  5,  this  angle  is  positive,  and  hence  the  first 
correction  applied  corresponds  to  a proper  (counterclockwise)  rotation 
in  the  missile  frame.  Were  Y|I  negative,  the  value  of  oq  would  be 
negative  and  a clockwise  rotation  in  the  missile  frame  would  be  indicated 
to  turn  the  missile  line  of  flight  toward  the  decoy. 

Having  determined  oq,  the  angle  8'  is  changed  according  to 

^new"  (12) 

Once  B/iew  "is  determined  using  (12),  the  line  of  flight  of  the  missile 
during  the  subsequent  time  increment  5t  can  be  calculated.  This  is  done 
in  effect  by  computing  the  equation  of  the  line  of  flight  from  the  current 
position  of  the  missile  (Xt,Y^)  to  the  revised  aim  point.  Referring  to 


Figure  6,  the  general  equation  for  this  line 
(in  the  horizontal  plane)  is 


(x^.y,.) 


Y = mX  + b 

The  slope  m is  given  by 


m * tan  (8*  + y’)  (14) 


I 


’‘s  ’‘int 


FIGURE  6 


[ while  the  y - axis  intercept  point  b is  determined  by  solving  (13)  with 

Y » Yt»  X = and  using  (14): 
b * Yt  - m Xt 

The  X - intercept,  which  is  the  revised  aim  point  of  the  missile,  is 
then  found  by  solving  (13)  for  Y=0: 

*t  (6new*f) 

The  missile  is  .then  advanced  along  the  new  line  of  flight,  and  the  ship 
and  chaff  cloud  positions  are  incremented.  If  the  missile  is  within  one 
time  increment's  distance  of  either  the  ship  or  the  decoy,  it  is  then 
determined  whether  the  missile  flight  path  will  take  it  over  the  ship. 
Otherwise,  the  entire  process  of  calculating  the  error  angle  oq  and 
moving  the  missile,  ship  and  chaff  cloud  is  repeated.  In  these  sub- 
sequent interactions  of  the  tracking  loop,  the  second  of  equations  (11) 
is  modified  by  substituting  X-j^t  (the  current  aim  point)  for  Xs. 

The  implementation  of  the  above  algorithm  in  computer 
software  requires  accounting  for  the  algorithm  for  computing  arctangents 
in  FORTRAN  IV  and  also  the  formulas  for  determining  the  rotation  angle 
of  the  missile  frame  and  the  aim  point.  Consequently,  it  was  necessary 
to  include  three  tracking  loops  to  deal  with  first  and  second,  third, 
and  fourth  quadrant  missile  locations. 
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IV.  EVALUATION  CRITERIA 

Implementation  of  the  models  in:  software  permits  calculation 
of  the  trajectories  of  the  ship,  decoy  and  missile  under  a wide  range  of 
initial  conditions.  This  in  turn  provides  the  information  required  to 
determine  the  relative  positions  of  the  objects  at  any  time,  and  thus  to 
evaluate  decoy  placement  effectiveness.  This  section  describes  the 
methodology  for  that  evaluation. 

The  simulation  program  evaluates  two  aspects  of  decoy  placement: 
whether  the  decoy  is  within  the  field  of  view  of  the  missile  for  a 
sufficient  time  to  affect  it  (given  that  the  missile  is  aimed  at  the 
ship),  and  assuming  that  this  is  so,  whether  the  decoy  will  cause  a 
sufficient  error  in  the  missile  to  cause  it  to  miss  the  ship.  It  is 
emphasized  from  the  start  that  the  second  of  these  criteria  cannot  be  as 
reliable  as  the  first,  because  it  is  based  on  the  assumption  that  track 
is  shifted  to  the  decoy  (and  not,  for  example,  to  the  ship-decoy  centroid) 
at  some  specific  time.  Nonetheless,  the  second  placement  criterion 
provides  an  indication  of  the  efficacy  of  the  decoy  in  causing  angle 
tracking  errors. 

The  first  success  criterion,  called  SUCCESS  1 (S]  for  brevity) 
is  based  on  the  purely  geometric  condition  that  the  decoy  is  within  the 
missile  field  of  view  for  a period  of  time  equal  to  or  greater  than  the 
decision  time  of  the  missile.  In  the  simulation,  time  is  stepped  in 
fixed  increments  of  5t  which  may  or  may  not  be  equal  to  the  decision 
time  To  at  the  option  of  the  user.  To  relate  fit  to  Tq,  a parameter 
called  ITS  is  defined  as  the  smallest  integer  satisfying  the  following 
relationship: 
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Then,  for  example,  if  Tq=  6t  ITS=2,  so  that  the  decoy  must  be  in  the 
field  of  view  at  two  successive  steps  in  time  to  count  as  a success  with 
respect  to  S].  The  implementation  of  this  condition  is  shown  in 
Figure  7,  which  indicates  the  structure  for  determining  the  outcome  of 
a single  engagement  in  the  simulation. 


FIGURE  7 

Determination  of  decoy  presence  within  the  missile  field  of 
view  (FOV)  is  accomplished  by  evaluating  two  quantities  referenced  to 
the  missile;  the  angle  from  the  ship  midpoint  to  the  most  distant 


cross-range  extremity  of  the  decoy,  and  the  extreme  distance  from  the 


ship  midpoint  to  the  decoy  boundary  along  the  missile  line  of  flight. 


If  the  ship  coordinates  in  the  missile  frame  are  (X5",Y5",Z5")  while 
those  of  the  decoy  are  (Xq",  Yq",  Zq"),  then  the  angle  CCANG  is  found 
be  solving  the  triangle  of  Figure  8 for  005  Y" 


using  the  law  of  cosines: 


CCANG  =cos 


D^-A^N  (15) 


where 


(x|:.Y».z||) 


S*  = X''2  + Y"2  + 7"2 

s s s 


= X"2  + y"2  + 7ti2 

C C 

A'  = CX''  -xp»  . CY"  + C2"  -z"y 

* c SC 


FIGURE  8 


To  the  angle  given  by  equation  (12)  is  added  the  angle  subtended  by 


half  of  the  decoy  as  viewed  by  the  missile.  If  DCC  is  the  radius  of  the 


chaff  cloud,  this  angle  is 

DANG  = tan  (DCC/D) 


That  is. 


ANGLE  = CCANG  + DANG 


ANGLE  < 6WHA 


where  BWHA  is  the  angle  from  missile  bores ight  to  its  3dB  one-way  antenna 
radiation  pattern,  then  the  decoy  is  considered  to  be  within  the  missile 


acceptance  angle. 


The  range  difference  between  the  ship  midpoint  and  the  decoy 


center  as  viewed  by  the  missile  is 


PRCC  = I S - D| 


To  this  we  add  DCC,  so  that  the  extreme  range  difference  between  the 


midpoint  of  the  ship  and  the  decoy  is 


RANGE  » PRCC  + DCC 


J 


Now  if. 


RANGE  < RGHW  (17) 

where  RGHW  Is  the  half-width  of  the  missile  range  gate,  then  the  decoy 
is  within  the  missile  range  gate.  If  (16)  and  (17)  hold  simultaneously, 
then  the  decoy  is  said  to  be  in  the  field  of  view  (FOV).  The  decoy  is 
successfully  placed  (S-|  is  true)  if  equations  (16)  and  (17)  hold 
simultaneously  for  ITS  successive  evaluations  subsequent  to  the  decoy's 
achievement  of  viability  (time  t>t[j+tf). 

Given  that  the  decoy  satisfies  $i,  the  missile  is  assumed  to 
immediately  initiate  maneuver  toward  the  decoy.  If  the  flight  of  the 
missile  toward  the  decoy,  projected  in  the  direction  of  flight  immediately 
prior  to  decoy  impact  or  as  the  missile  passes  the  ship,  crosses  any  part 
of  the  ship,  it  is  assumed  that  the  missile  will  hit  the  ship.  This  is 
evaluated  by  checking  the  signs  of  the  Y"  coordinates  of  the  ship's  bow 
and  stern.  If  they  are  both  of  the  same  sign,  i.e. 

Y^g>0andY^S>0 


or 

Y^g<0andY^5<0 


(18) 


then  the  second  success  criterion,  SUCCESS  2(or  Sg)  is  satisfied. 

In  summary:  $2  is  true  provided  that  equation  (18)  holds  at 
the  time  the  missile  impacts  the  decoy  or  passes  the  ship. 
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V. 


SIMULATION  SOFTWARE 


The  Decoy  Placement  Simulation  RFDPS  was  written  in  the  version 
of  the  FORTRAN  IV  (6)  language  utilized  by  the  DEC  System  10  computer 
(DEC,  1972).  The  program  was  designed  to  facilitate  its  use  with  time- 
shared  computer  facilities,  in  particular  those  of  TYMSHARE,  Inc. 

Although  FORTRAN  IV  (6)  is  a standard  language  and  time  sharing  computing 
procedures  and  capabilities  are  fairly  uniform,  it  is  impossible  to 
guarantee  that  the  program  described  in  this  section  can  be  utilized  on 
other  systems  without  alteration. 

In  order  to  execute  RFDPS,  it  is  necessary  to  compile  the 
program  and  create  two  input  data  files.  The  methods  for  accomplishing 
this  using  TYMSHARE  facilities  are  described  in  detail  in  XEXEC 
(TYMSHARE,  1974)  and  EDITOR  (TYMSHARE,  1969),  so  this  section  will 
concentrate  on  the  file  contents  rather  than  the  method  for  their  creation. 

The  first  of  the  input  data  files  is  referred  to  as  the 
parameter  file.  It  contains  the  parameters  listed  in  Table  1;  the 
layout  is  as  shown  in  Table  2.  In  accordance  with  good  programming  style 
(Kernighan  and  Plauger,  1974)  all  parameter  names  consist  entirely  of 
alphabetical  characters. 

As  indicated  in  Table  1,  three  seed  numbers  are  required  as 
input  parameters;  they  are  used  in  generating  the  true  wind  direction 
and  speed,  and  the  missile  initial  direction.  Three  seeds  are  required 
because,  as  shown  by  MacLaren  and  Marsaglia  (1965)  the  multiplicative 
congruential  generator  cannot  produce  uniformly  distributed  triples 
(See  also  McQuay  (1973)).  Since  the  seeds  and  the  multiplier  should  be 
relatively  prime  with  respect  to  one  another,  all  three  seed  numbers 
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NAME 

EXPLANATION 

RESTRICTIONS 

IRPT 

NUMBER  OF  CASES  TO  BE  INCLUDED  IN 
DETAILED  REPORT 

INTEGER  < IMMX 

1C 

VERSION  1C  = 0 GIVES  MONTE  CARLO  RUN 
1C  ^ 0 GIVES  DETERMINISTIC  RUN 

INTEGER 

IMMX 

NUMBER  OF  CASES  IN  RUN  (1C  = 0) 

NUMBER  OF  THREAT  BEARINGS  PER  VALUE  OF 
VW  AND  THETA  (1C  #0) 

INTEGER  <9999 

ISEEO  (1) 

SEEDS  FOR  RANDOM  NUMBER  GENERATORS 
FOR  VW,  THETA  AND  GAMMA,  RESPECTIVELY 

PRIME,  OF  FORM  8N  ± 3 
FOR  SOME  INTEGER  N 

WVL 

MAXIMUM  WIND  SPEED,  F/S  (MONTE  CARLO  VERSION) 

INTEGER  < 100 

DELTA! 

TIME  INCREMENT  BETWEEN  SUCCESSIVE 
EVALUATIONS  OF  PLACEMENT,  SECONDS 

REAL>0 

GAMZRO 

MINIMUM  VALUE  OF  THREAT  DIRECTION,  DEG 

0 < GAMZRO  < GAMMAX 

GAMMAX 

MAXIMUM  VALUE  OF  THREAT  DIRECTION,  DEG 

GAMZRO  < GAMMAX  < 360° 

TF 

DECOY  TIME  OF  FLIGHT,  S 

REAL>0 

TB 

DECOY  TIME  OF  SIGNATURE  DEVELOPMENT,  S 

REAL  >0 

BSHO 

HALF  LENGTH  OF  SHIP,  F 

REAL>0 

RCH 

HORIZONTAL  RANGE  OF  DECOY,  F 

REAL>0 

DCC 

RADIUS  OF  DECOY,  F 

REAL>0 

HZERO 

ALTITUDE  OF  DECOY  AT  TIME  OF  BURST,  F 

REAL>0 

VCV 

FALL  RATE  OF  DECOY,  F/S 

REAL 

DISTL 

LOCATION  OF  DECOY  REFERENCED  TO  AMIDSHIPS,  F 

REAL 

VSKN 

SHIP  SPEED,  KNOTS 

INTEGER  >0 

IPHI 

LAUNCHER  TRAIN  ANGLE  REFERENCED  TO  SHIP 
LONGITUDINAL  AXIS 

0 < INTEGER  < 180° 

RZERO 

MISSILE  RANGE  AT  TIME  OF  DECOY  LAUNCH,  F 

REAL>0 

ZTZERO 

MISSILE  ALTITUDE,  F 

REAL>0 

VT 

MISSILE  SPEED,  F/S 

REAL>0 

YKAPA 

MISSILE  NAVIGATIONAL  CONSTANT,  S'^ 

REAL>0 

BWHA 

MISSILE  3 dB  HALF  BEAMWIDTH 

REAL>0 

RGHW 

MISSILE  RANGE  GATE  HALF  WIDTH,  NS 

REAL>0 

ITS 

NUMBER  OF  TIMES  DECOY  MUST  BE  IN  FOV 

INTEGER  >0 

VWZRO 

FOR  IC¥=0,  THE  MINIMUM  WIND  SPEED,  KN 

REAL  >0 

VWI 

FOR  IC^O,  THE  WIND  SPEED  INCREMENT,  KN 

POSITIVE  REAL  SUCH  THAT 
VWM  = VWZRO  + KVWI  FOR 
SOME  INTEGER  K 

VWM 

FOR  IC^O,  THE  MAXIMUM  WIND  SPEED,  KN 

REAL>  VWO 

THZRO 

FOR  IC?tO,THE  MINIMUM  TRUE  WIND  DIRECTION,  DEG 

REAL,  0 < THZRO  < 360° 

THETI 

FOR  IC^O,  THE  TRUE  WIND  DIRECTION 
INCREMENT,  DEG 

POSITIVE  REAL  SUCH  THAT 
THETM  + THZRO  + NTHETI 
FOR  SOME  INTEGER  n 

THETM 

FOR  IC=)fcO,  THE  MAXIMUM  TRUE  WIND 
DIRECTION,  DEG 

REAL  >THETO 

I 


TABLE  1 


should  be  chosen  to  be  prime.  Further,  since  the  period  of  the 
generator  is  maximized  (and  equal  to  where  b is  the  word  length  of 
the  computer)  for  ISEED=8N+3  where  N is  an  integer  (Newman  and  Odell, 
1971),  the  seed  numbers  should  all  be  of  that  form.  Such  is  always 
possible,  and  can  be  accomplished  with  little  effort  (Carmichael,  1914). 
A partial  table  of  good  seed  numbers,  drawn  from  Abramowitz  and  Stegun 
(1973),  is  included  as  Appendix  B. 

The  second  of  the  input  files  contains  the  true  wind  speed 
distribution  table,  which  is  discussed  in  the  Appendix.  The  table 
DIST  contains  100  pairs  of  four  digit  integers  (the  program  will  accept 
150  pairs  without  modification).  The  first  entry  is  the  true  wind  in 
feet  per  second;  the  second  is  a four  digit  number  which  is  used  in 
conjunction  with  the  random  number  generator  for  the  true  wind  speed 
to  obtain  a sample  population  resembling  the  desired  parent  population. 

The  program  can  be  executed  using  DO  loops  to  generate  VW 
(true  wind  speed)  and  THETA  (true  wind  direction)  in  addition  to  the 
Monte  Carlo  method  of  generation.  This  mode  of  operation  can  be  of 
great  utility  when  a limited  range  of  wind  conditions  is  of  interest, 
as  when  evaluating  decoy  placement  against  a scenario  which  is  specific 
in  terms  of  geographic  or  weather  conditions  and  ship  direction  of 
advance. 

To  run  the  program  using  TYMSHARE,  the  following  sequence  of 
operations  is  required.  The  sequence  below  assumes  that  the  data  files 
have  been  written  and  the  program  compiled.  Underlined  statements  are 
entered  by  the  user;  the  symbol  t indicates  a carriage  return. 


(Log  on  to  system) 


- EXECUTE  RFDPS  <t  i 

i 

LOADING  i 

EXECUTION  I 

ENTER  PARAM  FILE  NAME:  PARAMF^  j 

j 

ENTER  NAME  OF  DIST  TABLE:  DIST(t 

- 

ENTER  SUM  RPT  FILE  NAME:  SUMARY(t 

ENTER  DTL  RPT  FILE  NAME:  DETAIL(t 

(Program  Executes)  | 

TOTAL  TRU  FOR  RUN=XXX  i 

i 

TRU  PER  CASE=YYY  j 

EXIT 

j 

-■ (Log  off  from  system  or  repeat  sequence  for  another  data  set) 

In  the  above  example,  the  file  PARAMF  contained  the  parameters  listed 

i 

in  lable  2,  DIST  was  the  table  of  true  wind  speed  probabilities;  the  | 

summary  file  SUMARY  and  the  detail  (case  by  case)  file  DETAIL  were  written. 

The  last  two  lines  gave  the  CPU  use  for  the  run  and  the  CPU  use  per  case.  I 

If  desired,  the  output  files  can  be  typed  on  the  remote 
terminal  simply  by  giving  the  command  TYPE  (File  Name)  t.  Alternatively, 
they  can  be  listed  (at  lower  cost)  on  a remote  high  speed  printer  by 
means  of  the  SPOOL  command.  Both  files  are  width  compatible  with  TTY 
size  paper  and  half  width  terminals. 

The  remaining  pages  of  this  section  contain  complete  listings  I 

j 

for  RFDPS  and  DIST.  The  main  program  requires  approximately  15.3k  bytes  | 

of  core  in  the  compiled  version  DIST  requires  an  additional  0.9k  bytes. 
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FIRST  LINE:  FORMAT  (71.  3F)  "PROGRAM  CONTROL" 

IRPT,  IMMX,  1C,  ISEED  (1),  ISEED  (2),  ISEED  (3),  VWL,  DELTAT, 
GAMZRO.GAMMAX 

SECOND  LINE:  FORMAT  (8F,  21)  “DECOY  SYSTEM/INSTALLATION 
AND  SHIP  PARAMETERS" 

TF,  TB,  BSHD,  RCH,  DCC,  HZERO,  VCV,  DISTL,  VSKN,  IPHI 

THIRD  LINE:  FORMAT  (6F,  I)  “MISSILE  PARAMETERS" 

RZERO,  ZTZERTO,  VT,  YKAPA,  BWHA,  RGHW,  ITS 

FOURTH  LINE:  FORMAT  (10A4)  "FILE  NAME" 

FILE  NAME 

FIFTH  LINE:  FORMAT  (6F)  "WIND  PARAMETERS" 

VWZRO,  VWI,  VWM,  THETZRO,  THETI,  THETM 


RFDPS  LISTING 


DOUBLE  PRECISION  NAME , DNAME , SNAME , FNAME 
INTEGER  SUCS1,SUCS2,W/L,VSKN 

DIMENSION  IVr.vT(2,150)  , ISEED  (3)  , IG  (24)  , IGF  (2,24)  , 
-GS(2,24)  ,IRW(24)  , IRV.T  (2,24)  ,RWS(2,24)  ,IW(13)  , 

-IVWF (2,13) ,VWS(2,13) ,IGLO(24) ,IGHI(24) , 

-IVLO(12) ,IVHI (12) ,IDLO(24) ,IDhI(24) ,IDSC(10) 
XSTRT=TRU(XXX) 

C***  READ  IN  PAI?AM  CARDS 
TYPE  4 

4 FORMAT  ('  ENTER  PARAM  FILE  NAI>IE:  ',$) 

ACCEPT  7,FNAI^J: 

7 FORMAT (AlO) 

OPEN  (8  , FNAI^Z , INPUT , SYMBOLIC , ERR=5) 

GO  TO  10 

5 TYPE  6 

6 FORMAT  (•  NO  PARAI4  FILE — PROGRAM  ABORTED.') 

STOP 

10  READ  (3,11)  IRPT,m IX, IC, ISEED, W.Ti,DELTAT,GAr'ZP.O,GAM?’i\X 

11  FOPJ'IAT(7I,3F) 

READ (8,12)  TF , TB , BSHD , RCH , DCC , HZERO , VCV , D ISTL , VSKN , IPKI 

12  FORMAT (8F, 21) 

READ  (8,13)  RZERO,  ZTZERO  , VT  , YKAPA  , BWUA  , RGIIW  , ITS 

13  FORflAT(6F,I) 

READ(8,14)  (IDSC(I) ,1=1,10) 

14  FORMAT (10A4) 

READ  (8,15)  WZ  RO , W7I , VWM , THZ  RO , THET I , THETM 

15  FORMAT (6F) 

CLOSE (8) 

C***  READ  IN  VW  DISTRIBUTION  TABLE. 

1=1 

TYPE  18 

18  FOPJIATC  enter  FILE  NATIE  OF  DIST  TABLE;  ',$) 

ACCEPT  7, NAME 

OPEN ( 9 , NAME , INPUT , SYMBOLIC , ERR=2 1 ) 

GO  TO  19 
21  TYPE  17 

17  FORI'IATC  no  W distribution  table  - PROG.  ABORTED.') 
STOP 

19  READ  (9,16,END=20)  IVWT (1 , I) , IVWT (2, I) 

16  FORMAT (214) 

1=1+1 

IF(I.GT.150)  GO  TO  20 
GO  TO  19 

20  IVTr-LX=I-l 
CLOSE (9) 

DO  25  J=1,IVTMX 
IF(VWL.LE.I\A\T(1,J)  ) GO  TO  30 
25  CONTINUE 
30  IM70C=IVWT’(2,J) 

C***  INITILIZE  COUNTERS 


ITHM=  (THETM-TIIZRO)  /THETI+1 . 0 
IVWM=  (VW1-1-W2RO)  /WI+1 . 0 

c***  ope::  i/o  files. 

TYPE  37 

37  FORMAT C*  ENTER  SUM  RPT  FILE  NAME:  ',$) 

ACCEPT  7,SNAI’IE 

OPEN  (10,  SNAIIE , OUTPUT , SYMBOLIC ) 

IF(IRPT.GT.O)  TYPE  38 

38  FORMAT  (•  ENTER  DTL  PJ>T  FILE  NAME:  ',$) 
IF(IRPT.GT.O)  ACCEPT  7,DNAIffi 

IF  (I RPT . GT . 0 ) OPEN  ( 11 , DNAME , OUTPUT , S YI-30LIC ) 

CRD=57. 29578 

CNMF=6076 

VS=VSKN*  CNMP/3  600 

PHI=IPHI/CRD 

SPHI=SIN(PHI) 

CPUI=C0S(PH1) 

BWHA=BWHA/CRD 
RGIIV;=RGIM*.9b4 
YDIS=VT*DELTAT 
C*****CLEAR  ARRAYS 
DO  40  1=1,12 
IG(I)=0 
IWV(I)=0 
40  IWV(13)=0 
DO  45  1=1,24 
IRW(I)=0 
45  IFAIL1=0 
IFAIL2=0 
IFAIL=0 

C*****LOOP  THRU  IM 

50  DO  2375  IM=1,IMMX 

Qlt***  ******  ft********* 

C IF  DETERMINISTIC,  LOOP  THROUGH  THETA  AND  VW 

Qit***  *******  ********* 

IF(IC.EQ.O)  GO  TO  60 
DO  2350  ITHET=1,ITHM 
THETA=THZRO+ (ITHET-1) *THETI 
THETA=THETA/CRD 
DO  2300  IVT7=l,IVtVM 
VW-VWZRO+  (IVW-1)  *VWI 
VW=VW*3600/CNMF 
STHETA=SIN (THETA) 

CTHETA=COS (THETA) 

60  CONTINUE 

76  IF(VWL.GT.O.O)  GO  TO  80 
GO  TO  110 

C*****GET  RANDOM  VW 

IF(IC.GT.O.O)  GO  TO  120 
80  YFL=RAN(ISEED(D) 


ITEMP=YFL* IKAX+ . 5 
IF (ITEMP.lt. 1)  GO  TO  80 
85  DO  90  I=l,IVTt!X 

IF(ITEMP.GT.IW?T(2,I) ) GO  TO  90 
GO  TO  95 
90  CONTINUE 

95  IF(IVIVT(1,I)  .GT.VV/L)  GO  TO  80 
VW=IWJT(1,I) 

ivw=iwrr(i,i) 

IF(VW.EQ.O.O)  GO  TO  110 
C***  GET  RANDOM  THETA. 

100  YFL=RAtJ(ISEED(2) ) 

K-YFL*379.0 

IF(K.LT.10.OR.K.GT.369)  GO  TO  100 

ITHETA=K-10 

THETA=ITHETA/CRD 

STHETA=SIN (THETA) 

CTHETA=COS (THETA) 

GO  TO  120 
110  ITHETA=0 
THETA==0 . 0 
STHETA=0.0 
CTHETA=1.0 
VW*0.0 

C***  GET  RANDOM  GAMI-IA. 

120  YFL=RAN(ISEED(3) ) 

K=YFL*  (GAIi^lAX+19.  0) 

IF(K.LT. (GAMZRO+10.0) .OR.K.GT. (GAMMAX+10.0) ) GO  TO  120 

IGAMMA=K-10 

GAMHA=  IG  AIIMA/C  RD 

SGAMMA=SIN  (GAMJ-IA) 

CGAfU4A=COS  (GAIIMA) 

C*****CALCULAT£  TEND 

BETA=ASIN  (VS*  SGA^m/VT) 

SBG»SIN (BETA+GAMMA) 

CBG=COS  (BETA+GAim) 

IF(IGAMMA.NE.0.AND.IGAMI4A.NE.180)  GO  TO  190 
IF  CIG/^JIA . EQ . 1 8 0 ) TEND=R2ERO/  (VT-VS ) 

IF ( IGAMMA . EQ . 0 ) TEND=R2ERO/ (VT+VS ) 

GO  TO  191 

190  TEND=RZERO*SGA^U4A/(VT*SLG) 

C***  CALCULATE  RWA 
191  IRWA=0.0 

IF(VWL.EQ.O)  GO  TO  250 
200  TEMP=VW*CTHETA+VS 

IF(ITHETA.EQ.i80)  GO  TO  210 
IF  (TEMP. NE. 0.0)  GO  TO  220 
IF(ITHETA.LT.180)  GO  TO  205 
ZETA=270.0 
GO  TO  240 


205  ZETA=yO.O 
GO  TO  240 

210  IF  (TE14P.lt.  0.0)  GO  TO  215 
ZETA=0. 0 
GO  TO  240 
215  ZETA=180.0 
GO  TO  240 

220  TE14P=  ( VW*  STHETA ) /TEMP 

IF  (TE:1P  . LT  . 0 . 0 ) TEMP=-TEMP 
ZETA=ATAN (TEMP) 

ZETA=ZETA*CRD 
TEMP=VW*CTKETA 
IF (TEMP . LT . 0 . 0 ) TEMP=-TEMP 
IF(ITHETA.GE.0.AND.ITHETA.LE.90)  GO  TO  240 
IF(ITHETA.LE.1&0)  GO  TO  230 
IF(ITHETA.LE.270)  GO  TO  235 
225  ZETA=360-ZETA 
GO  TO  240 

230  IF (TEMP. LE. VS)  GO  TO  240 
ZETA=180-ZETA 
GO  TO  240 

235  IF (TEMP. LE. VS)  GO  TO  225 
ZETA=ZETA+180 
240  IRWA=ZETA 
C***  CALCULATE  VRW 

250  VRW=SQRT  CVW*W';+VS*VS+2*VW*VS*CTHETA) 

C***  INITILIZE  TIME. 

300  T=TB+TF 

P1=RCH*CPHI+DISTL+VS*TF 

P2=VW*CTHETA 

P3=RCH*SPKI 

P4=VW* STHETA 

P5»RZERO*SGAMMA 

P6=RZERO*CGAKI>lA 

ICNTS=0 

sues 1=0 

sues 2=0 

TSUCS1=0 

C**** ‘CALCULATE  COORDS  IN  EARTH  FRAME 
320  XC=P1-T*P2 
YC=P3-T*P4 
ZC=HZERO-VCV* (T-TF) 

IP(ZC.LT.O.O)  ZC=0.0 
XS»VS*T 

XT«P6+VS*T-T*P6/TEND 
326  yr=P5* (1-T/TEND) 

C***‘*CALCULATE  COORDS  IN  TRACKER  FRAME 
XDS=  (XT-XS) *CBG+YT*SLG 
YDS=  (XS-XT)  *S13G+YT‘CBG 
XDSB- (XT-XS-BSHD) *CBG+YT*SBG 
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XDSS=  (XT-XS+BSIID)  *CBG+YT*SBG 
YDSB= (XS+BSHD-XT) *SBG+YT*CBG 
YDSS= (XS-BSHD-XT) *SBG+YT*CBG 
RHO=SQRT (XDS*  *2+YDS**2 ) 

ZT=ZTZERO 
328  ZDS=-ZT 

XDC= (XT-XC) *CBG+ (YT-YC) *SBG 
YDC= (XC-XT) *SBG+  CYT-YC) *CBG 
ZDC=ZC-ZT 

C***  CALCULATE  CCAl^G. 

600  AH=XDS**2+YDS**2+ZDS**2 

BH=XDC  *XDC+YDC*YDC+ZDC*ZDC 

CH= (XDS-XDC)* (XDS-XDC)  + (YDS-YDC)* (YDS-YDC)  + (ZDS-ZDC) **2 
Z3= (AH+BH-CK) / (2*SQRT {AH*BH) ) 

IF(Z3,GE.1.0)  Z3=1.0 
CCANG=ACOS(Z3) 

C*****CALCULATE  RANGE  DIFFERENCE 
800  PRCC=ABS (SQRT (AH) -SQRT (BH) ) 

C*****EVALUATE  SUCCESS  ONE 
850  DANG=ATAN (DCC/SQRT (EH) ) 

ANGLE=CCANG+DANG 

RANGE=PRCC+DCC 

IF  (ANGLE. LE.BWEA. AND. RANGE. LE.RGm7)  ICNTS=ICNTS+1 
IF(ICNTS.GE.ITS)  SUCSl=l 
IF(SUCSl.EQ.l)  GO  TO  890 
810  T=T+DELTAT 

IF(T.LE.TEND)  GO  TO  320 
GO  TO  1000 

C***  TRACK  DECOY  AND  EVALUATE  SUCS2 
890  TSUCS1=T 

IF(XT.LT. 0.0. AND.YT.lt. 0.0)  GO  TO  930 
IF(XT.GT.O.O.AZJD.YT.LT.O.O)  GO  TO  960 
895  GAI1PR=ATAN(YT/XT) 

IF(GAMPR.LT.O.O)  GAMPR=GANPR+180/CRD 
BETPR=ASIN(XS*SIN(GAIIPR)/SQRT(  (XT-XS)**2+YT**2) ) 

900  SGBP=SIN(G.A>'PR+BETPR) 

CGBP=COS (GAMPR+LETPR) 

XDS= (XT-XS) *CGBP+YT*SGBP 
XDC= (XT-XC) *CGBP+ (YT-YC) *SGBP 
YDC= (XC-XT) *SGBP+ (YT-YC) *CGBP 
IF(XDC.LT.YDIS.OR.XDS.LT.YDIS)  GO  TO  920 
ALPHD=ATAN (YDC/XDC) 

BETPR=BETPR+YKAPA*ALPED*DELTAT 
CGBP=COS  (GAIIPR+3ETPR) 

SGBP=SIN (GAMPR+BETPR) 

XINT=XT-YT*CGBP/SGBP 

XT=XT-YDIS*CGBP 

YT=YT-YDIS*SGBP 

XC=XC-P2*DELTAT 

YC*=YC-P4*DLLTAT 
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910  XS=XS+VS*DELTAT 
GAtlPR=ATAII  (YT/XT) 

IF  (GAMPR.lt.  0.0)  GAI1PR=GAMPR+13 0/CRD 

BETPR=ASIN  (XIWT*SIN  (GAflPR)  /SORT  ( (XT-XINT)  **2+YT**2) ) 

GO  TO  900 

920  YDSB=  (XS+BSIiD-XT)  *SGBP+YT*CGBP 
YDSS=YDSli-2*BSiiD*SGBP 

IP (YDSB. GT. 0. 0. AND. YDSS.GT. 0.0. OR.YDSB.lt. 0.0. AND. 
-YDSS.LT.0.0)  GO  TO  990 
GO  TO  900 

C*****TRACKING  LOOP;  3RD  QUADRANT 
930  GAI'IPR=ATAN  (YT/XT) 

BETPR=ASIN  (XS*SIK  (GAIIPR)  /SQRT  ( (XT-XS)  **2+YT**2)  ) 

935  SGBP=SIN(GAtlPR“BETPR) 

CGBP=C0S  (GAIIPR-BETPR) 

XDC= (XC-XT) *CGBP+ (YC-YT) *SGBP 
YDC= (XT-XC) *SGBP+ (YC-YT) *CGBP 
XDS= (XS-XT) *CGEP-YT*SGBP 

IP(XDC.LT.YDIS.OR.XDS.LT.YDIS.OR.YT.GT.-YDIS)  GO  TO  940 
ALPHD=7iTAN  (YDC/XDC) 

BETPR=BETPR-YKAPA*ALPHD*DELTAT 
CGBP=COS (GAMPR-BETPR) 

SGBP=SIN  (GAIIPR-BETPR) 

XINT=XT-YT*CGBP/SGBP 

XT=XT+YDIS*CGBP 

YT=YT+YDIS*SGBP 

XC=XC-P2*DELTAT 

YC=YC-P4*DELTAT 

XS=XS+VS*DELTAT 

GAMPR=ATAN (YT/XT) 

BETPR=ASIN  (XINT*SIN (GAMPR)/SQRT ( (XT-XINT) **2+YT**2) ) 
IP(XT.LT. 0.0. AND.YT.lt. 0.0)  GO  TO  935 
IP(XT.GT. 0.0. AND.YT.lt. 0.0)  GO  TO  960 
GO  TO  095 

940  YDSB= (XT-XS-BSHD) *SGBP+YT*CGBP 
YDSS=YDSB+2*BSUD*SGBP 

IP (YDSB. GT. 0.0. AND. YDSS.GT. 0.0. OR.YDSB.lt. 0.0. AND. 
-YDSS.LT.0.0)  GO  TO  990 
GO  TO  980 

C*****TRACKING  LOOP:  4TH  QUADRANT 
960  GAMP R=-ATAN (YT/XT) 

BETPR=ASIN(XS*SIN(GAMPR) /SQRT ((XT-XS) **2+YT**2) ) 

965  SGBPsSIN  (GAMPR-t-BETPR) 

CGBP-COS (GAMPR+BETPR) 

XDC» (XT-XC) *CGBP+ (YC-YT) *SGBP 
YDC« (XT-XC) *SGbP+ (YT-YC) *CGBP 
XDS« (XT-XS) *CGBF-YT*SGBP 

IF(XDC.LT.YDIS.OR.XDS.LT.YDIS.OR.YT.GT.-YDIS)  GO  TO  970 
ALPHD-ATAN (YDC/XDC) 

BETPR-BETPR-YKAPA*ALPHD*DELTAT 
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CGBP=COS  (GAIlPR+iiETPR) 

SGLP=SIN  (C7J1PR+BETPR) 

XINT=XT+YT*CGEP/SGBP 

XT=XT-YDIS*CGBP 

YT=YT+YDIS*SGBP 

XC=XC-P2*DELTAT 

YC=YC-P4*DLLTAT 

XS=XS+VS*DELTAT 

GAt-IPR=-ATAN  (YT/XT) 

BETPR=ASIN (XINT*SIN (GAMPR) /SQRT ( (XT-XINT) **2+YT**2) ) 
IF(XT.GT.O.O.AI'ID.YT.LT.O.O)  GO  TO  96b 
IF(XT.LT.0.0.7iliD.YT.LT.0.0)  GO  TO  930 
GO  TO  895 

970  YDSB= (XT-XS-BSHD) *SGBP+YT*CGBP 
YDSS=YDSB+2*BSIID*SGBP. 

IF  (YDSB.GT.O.O.AIID.YDSS.GT.  0.0. OR.YDSB.lt.  0.0. MID. 

“YDSS.LT.O. 0)  GO  TO  990 
980  SUCS2=0 

GO  TO  1000 
990  SUCS2=1 

C***  WRITE  DETAIL  REPORT 
1000  CONTIIIUE 

IG7^.MA=GM ’J-iA*  CPD 
IVRW=VRW*  360 0/C!II  IF+ . 5 
IF(IM.GT.IRPT)  GO  TO  2010 
IFCIM.GE.2)  GO  TO  1100 
WRITE  (11,1050)  FMM1E,VSKN,IPHI 
1050  FORMAT (IX, A10,2X, 'DETAILED  REPORT  FOR  VS  = ',12, 

-'KNOTS  M'lD  PHI  = ',  13 ,' DEGREES '/ 

-4X,  'GMUIA'  ,5X,  'RWA'  ,5X,  'VRW'  ,2X,  'SUCSl'  ,2X,  'TSUCSl'  ,2X,  'SUCS2'/ 
-3X, 'DEG.REL.  DEG.REL.  KlI*  ,10X, 'SEC. '//) 

1100  WRITE  (11,1150)  IGATH  ‘A , I RWA , I VRW , SUCS 1 , TSUCS 1 , S UCS  2 
1150  FORiyiAT(5X,I3,6X,I3,5X,IJ,3X,Il,5X,F4.1,5X,Il) 

C***  ACCUMULATE  RUN  STATISTICS 
2010  IF  (SUCSl. EQ.O)  IFAIL1=IFAIL1+1 

IF  (SUCSl . EQ . 1 . AND . SUCS2 . EQ.O)  IFAIL2=«IFAIL2+1 
IF  (SUCSl . EQ . 0 . OR . SUCS  1 . EQ . 1 . 7U4D  . SUCS 2 . EQ . 0 ) IFAIL=IFAIL+ 1 
C***  ACCUMULATE  THREAT  GEOMETRY  AND  REL.  WIND  7VNGLE  STATISTICS 
DO  2050  1=1,24 
M»15* (I-l) 

N-15*I 

IF(M.LE.IG7^1MA.AND.N.GT.IG7^'A)  IG  (I)  =IG  (I) +1 
IF  (M . LE . IG A1>1MA . AND . N . GT . IGAMMA . AND . SUCSl . EQ . 0 ) 
-IGF(1,I)=IGF(1,I)+1 

IF  (M . LE . IGAMMA . AND . N . GT . IG7U^I^ A.  AND . SUCS 2 . EQ . 0 . 

-AND. SUCSl. EQ.l)  IGF (2 , I) =IGF (2 , I) +1 
2050  CONTINUE 

DO  2100  1=1,24 
M»15*(I-1) 

N-15*I 
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IF  (M  . LE  . I P V7A . A WD . I'J . GT . I RWA ) I RW  ( I ) = I RW  ( I ) + 1 

if{m.le.irv;a.and.n.gt.irwa.aiid.sucsi.eq,o) 

-IRWF(1,I)=IRV7F(1,I)+1 
IF (M . LE . I RWA . AND . N . GT . IRWA . AND . SUCS  2 , EQ . 0 . 
-AND.SUCSl.EQ.l)  IRWF (2,I)=IRVrt' (2,1)+! 

2100  CONTINUE 

DO  2200  1=1,12 
M=5*(I-1) 

N=5*I 

IF (M . LE . I VRW . AND . N . GT . I VRW)  I W ( I ) = IWV ( I ) +1 
IF  (M , LE . IVRW . AND . N . GT . IVRW , AIJD . SUCSl . EQ . 0 ) 
-IVtVF  (1 , 1)  =IVV<T  (1 , 1) +1 
IF (M . LE . IVRW . AND . N . GT . IVRW . AND . SUCS  2 . EQ . 0 . 
-AND.SUCSl.EQ.l)  IVWF  (2 , 1)  =IV\7F  (2 , 1) +1 
2200  CONTINUE 

IF (IVRW. LT. 60)  GO  TO  2250 
IWV{13)=I\'TV{13)+1 

IF  (SUCSl. EQ.O)  IVt7F(l,13)=IV;7F(l,13)+l 
IF(SUCS2.EQ.0)  1V\‘1F  (2 , 13)  =IVWF  (2 , 13) +1 
2250  IF(IC.EQ.O.O)  GO  TO  2375 
C***  END  VI'/  LOOP 


2300  CONTINUE 
C***  END  THETA  LOOP 
2350  CONTINUE 
C***  END  IM  LOOP 
2375  CONTINUE 
C***  WRITE  SUi'lMARY  REPORT 
2400  DO  2450  1=1,24 

IF(IG(I)  .NE.O)  GO  TO  2425 

GS(1,I)=0 

GS(2,I)=0 

GO  TO  2450 

2425  GS(l,I)=FLOAT(IG(I)-IGF(l,I)  )/IG(I) 

IF(IG(I)  .EQ.IGF(1,I) ) GO  TO  2450 

GS (2 , I) =FLOAT (IG (I ) -IGF (2,1) -IGF (1 , I) ) / (IG (I) -IGF (1,1)) 

2450  CONTINUE 

DO  2500  1=1,24 

IF(IRWd)  .NE.O)  GO  TO  2475 

RWS(1,I)=0 

RWS(2,I)=0 

GO  TO  2500 

2475  RWS(l,I)=FLOAT(IRW(I)-IRl>7F(l,I)  )/IRW(I) 

IF(IRWd)  .EQ.IRWF(1,I)  ) GO  TO  2500 

RWS (2 , I) =FLOAT (IRW (I ) -IRWF (2 , I) -IROT (1,1))/ (IRW (I) -IRWF (1,1)) 
2500  CONTINUE 

DO  2600  1=1,13 

IF (IWV (I) .NE.O)  GO  TO  2550 

W/S(1,I)=0 

VWS(2,I)»0 

GO  TO  2600 
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2550  W/S(l,I)=FLOAT(HW(I}-IVWF(l,I)  )/IVA7(I) 

IF(I;TV(I)  .EQ.IVhTd,!) ) GO  TO  2600 

VWS  C2 , 1)  =FL0AT  (IV, V (I)  -IVV-T  (2 , 1)  “IW;F  (1,1))/  (IWV (I)  -IVWF  (1,1)) 
2600  CONTINUE 

SRl=FLOAT  (I^U•IX--IFAIL1)  /IMTIX 
IF(IFAILl.EQ.IMIIX)  GO  TO  2650 
SR2=FLOAT  (II41-iX-IFAIL2-IFAILl)  / (IMMX-IFAILl) 

2650  SR3=FLOAT  (IMHX-IFAIL1-IFAIL2) /II-IMX 
BSUD2=2*ESHD 
WRITE (10, 2700)  FNAME 
2700  FORMAT (15X, 'INPUT  FILE: ' ,2X,AlO) 

WRITE (10, 2710)  VSKN 

2710  FORMAT (15X, 'SHIP  SPEED: ', 4X,I2, IX, 'KNOTS ' ) 

WRITE(10,2720)  IPHI 

2720  FORMAT (2X, 'LAUNCHER  POINTING  ANGLE :', 3X, 13 ,1X, 'DEG  REL') 

WRITE (10, 2730)  ESHD2 

2730  FORMAT (14X, 'SKIP  LENGTH: ', 3X,F4 . 0 , IX, 'FEET ' ) 

WRITE (10, 2740)  DISTL 

2740  FORMAT  (8X, 'LAUKCilLR  LOCATION:  ' ,2X,F5.1,1X, 'FT  FM  CTR'/) 

WRITE  (10, 2750)  I!i;iX 

2750  FOR!'!AT(2X, 'RUN  STATISTICS:  TOTAL  NUMBER  OF  RUNS  : ',14) 

WRITE (10,2760)  IFAILi 
2760  FOR2N^T(19X, 'TOTAL  NUMBER  OF  FAILl 
WRITE (10, 2770)  SRi 

2770  FORMAT (19X, 'TYPE  1 SUCCESS  RATIO 
WRITE (10, 2780)  IFAIL2 
2780  FORi'IAT(19X, 'TOTAL  NUMBER  OF  FAIL2 
WRITE (10, 2790)  SR2 

2790  FOR!N'.T(19X, 'TYPE  2 SUCCESS  iLATIO 
WRITE (10, 2795)  SR3 

2795  FOPJIAT(19X, 'CENTROID  SUCCESS  RATIO 
WRITE (10, 2800) 

2800  FORMAT (6X, 'FAILURES  BY  TRACKER  DIRECTION'/ 

-2X, 'DIRECTION' ,2X, 'NO.OF' ,4X, 'FAILl* ,7X, 'FAIL2'/ 

-2X, • (DEG  REL) ' ,2X, 'CASES', 2X, 'NO. ',2X, 'FRAC. ',2X, 'NO. ' , 

-2X, 'FRAC. '//) 

DO  2820  1=1,24 
IGLO(I)=15* (I-l) 

IGIII(I)=15*I 

WRITE(10,2810)  IGLO(I)  ,IGHI(I)  , IG  (I)  , IGF  (1 , 1)  ,GS(1,1)  , 
-IGF(2,I)  ,GS(2,I) 

2810  FORMAT(3X,I3, ,I3,3X,I4,2X,I4,2X,F4.2,2X,I4,2X,F4.2) 

2820  CONTINUE 

WRITE (10, 2840) 

2840  FORMAT (//3X, 'FAILURES  BY  RELATIVE  WIND  DIRECTION'/ 

-2X, 'DIRECTION  NO.OF' ,4X, 'FAILl' ,7X, 'FAIL2'/ 

-2X,'(DEG  REL)  CASES  NO.  FRAC.  NO.  FRAC.'//) 

DO  2860  1=1,24 
IDLO  (I)»15*(I-1) 

IDHI  (I)=15*I 


',14) 

*,F4.2) 

',14) 

',F4.2) 

•,F4.2/) 
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WRITE  (10, 2a50)  IDLO(l)  ,IDUI  (I)  , IRW (I)  , IRViF  (i , I)  ,nV7S(l,I)  , 
-IRV.T(2,I)  ,KKS  (2,1) 

2850  F0RI1AT(JX,I3,  ,I3,3X,I4,2X,I4,2X,F4.2,2X,I4,2X,F4.2) 

2860  CONTINUE 

WRITE(10,2870) 

2870  F0R11AT(//4X, 'FAILURES  BY  RELATIVE  WIND  SPEED'/ 

-3X,  'SPEED  IIO.OF'  ,4X,  'FAILl'  ,7X,  'PAIL2'/ 

-3X,' (KN) *,3X, 'CASES  NO.  FPAC.  NO.  FRAC. '//) 

DO  2890  1=1,12 
IVLO(I)  = 5*  (I-l) 

IVUI (I)=5*I 

WRITE  (10,2380)  IVLO(I)  , IVHI  (I)  ,IIW(I)  ,IV\^F(1,I)  ,VWS(1,I)  , 
-IVWF(2,I)  ,\AVS  (2,1) 

2880  FORMAT(3X,I2, '-' ,I2,2X,I4,2X,I4,2X,F4.2,2X,I4,2X,F4.2) 

2890  CONTINUE 

WRITE  (10, 2900)  IVA/(13)  ,IWJF(1,13)  ,VWS(1,13)  , I\A'7F  (2 ,13 ) ,WvS(2,13) 
2900  FORMAT(3X, 'GE  60 ' , 2X, 14 , 2X, 14 , 2X,F4 . 2, 2X, 14 , 2X,F4 . 2/) 

C###### 

XNOW=TRU (XXX) -XSTRT 
TYPE  4006,XKOW 

4006  FORMAT ('  TOTAL  TRU  FOR  RUN  = ',F8,4) 

XNO=XNOW/IIUIX 

TYPE  4007,XKO 

4007  FORi^iATC  TRU  PER  CASE  = *,F10.7) 

C#«####### 

STOP 

END 
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APPENDIX  A 


THE  DISTRIBUTION  OF  TRUE  WIND  SPEEDS  AT  SEA  AND 
AN  ALGORITHM  FOR  CALCULATING  RELATIVE  WIND 
A.l  DISTRIBUTION  OF  TRUE  WIND  SPEEDS 

The  decoy  placement  simulations  described  in  this  volume  require 
input  data  concerning  true  wind  conditions  since  these,  along  with  the 
ship's  speed,  determine  the  relative  motion  of  a free  falling  airborne 
decoy  with  respect  to  the  ship,  and  hence  also  with  respect  to  the 
missile.  Whereas  the  direction  of  the  true  wind  can  be  assumed  to  be 
uniformly  distributed  to  good  approximation  in  modeling  ensembles  of 
encounters,  the  true  wind  speed  cannot.  Hence,  it  was  desirable  to 
obtain  empirical  data  on  true  wind  speed  frequencies  of  occurence  in  a 
maritime  environment. 

Fortunately,  such  data  were  reported  by  Long  and  his  co-workers 
(Long,  et  al,  1965).  Their  report  contained  a curve  giving  true  wind 
speeds  and  their  frequencies  of  occurrence;  the  curve  resulted  from  the 
aggregation  of  data  for  a large  number  of  locations  in  several  oceans 
over  the  period  of  many  years. 

Long's  curve  was  analyzed  to  determine  if  it  was  reasonably 
well  approximated  by  any  known  probability  density  function.  The  results 
of  the  analysis  indicated  an  excellent  fit  for  a slightly  modified 
distribution  having  eight  degrees  of  freedom.  The  words  "slightly 
modified"  are  significant  here,  even  though  the  modification  was  extremely 
slight*.  It  must  be  emphasized  that  no  claim  is  being  made  for  a 
distribution  with  8 degrees  of  freedom  being  the  result  of  some  underlying 
property  of  nature.  We  could  have  done  just  as  well  with  a polynomial 
curve  fit. 

* It  consisted  of  asigning  a non-zero  probability  (P=10"^)  to  VW=0 


Once  a satisfactory  fit  to  Long's  curve  was  obtained,  it  was 
quantized.  The  result,  plotted  again  in  the  form  of  a curve,  is  shown 
in  Figure  Al.  Arranged  in  tabular  form,  the  curve  yields  the  table 
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A2.  CALCULATION  OF  RELATIVE  WIND  SPEED  AND  ANGLE 


Given  the  direction  of  the  wind  0 as  measured  in  the  earth  frame 


the  wind  speed  and  the  ship  speed 
V5,  the  relative  wind  speed  and 
the  relative  wind  direction  3 can  be 
determined  by  solving  the  vector  dia- 
gram shown  in  Figure  A. 2.  Using  the  law 
of  cosines. 


Vrw=(V*S+V^W+2VwVsCOS6)’'*  A.l 

Using  the  law  of  sines, 

3=sin”^(Vwsine/VRvi)  A. 2 

Now  equation  (A.l)  uniquely  specifies  Vrw  for  any  triple 
<Vn.  Vs,  e);  however  (A. 2)  is  multivalued:  for  given  Vy,  Vrw  and 

9 /N+l)ir/2,  it  admits  roots  3l  and  32  such  that  32=tt-3i  . Clearly,  only 
one  of  the  roots  can  be  physically  meaningful,  but  ostensibly  it  could  be 
either  one.  For  example,  if  9=0,  3 will  be  equal  to  0°  if  Vs>Vy,  and 
equal  to  180°  if  Vy>V£.  The  generalization  of  this  observation,  along 
with  recognition  of  the  fact  that  3 will  always  be  in  the  hemicircle 
containing  0,  V5,and  -Vs  allows  design  of  the  algorithm  for  calculating 
the  relative  wind  angle.  The  rules  for  calculation  of  3 are  as  follows: 


1)  If  0 lies  in  the  first  quadrant,  so  does  B 

2)  If  e lies  in  the  second  quadrant,  B lies 

in  the  second  quadrant  if 

|VyCos0l>Vs; 

otherwise  it  lies  in  the  first  quadrant 

3)  If  0 lies  in  the  third  quadrant,  B is  also  j 

'in  the  third  if  - 

|Vwcos0|>V5;  " 

otherwise  B lies  in  the  fourth  quadrant  I 

4)  If  0 lies  in  the  fourth  quadrant,  B is  also 

in  the  fourth  quadrant. 

I 

i 

i 

f • ■ : 
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FIVE  DIGIT  PRIMES  OF  FORM  PRIME  = 8N+3 
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